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An acoustic analogy is developed to predict the noise from jet flows. It contains two 
source models that independently predict the noise from turbulence and shock wave shear 
layer interactions. The acoustic analogy is based on the Euler equations and separates 
the sources from propagation. Propagation effects are taken into account by calculating 
the vector Green’s function of the linearized Euler equations. The sources are modeled 
following the work of Tam and Auriault, Morris and Boluriaan, and Morris and Miller. 
A statistical model of the two-point cross-correlation of the velocity fluctuations is used 
to describe the turbulence. The acoustic analogy attempts to take into account the cor- 
rect scaling of the sources for a wide range of nozzle pressure and temperature ratios. 
It does not make assumptions regarding fine- or large-scale turbulent noise sources, self- 
or shear-noise, or convective amplification. The acoustic analogy is partially informed by 
three-dimensional steady Reynolds- Averaged Navier-Stokes solutions that include the noz- 
zle geometry. The predictions are compared with experiments of jets operating subsonically 
through supersonically and at unheated and heated temperatures. Predictions generally 
capture the scaling of both mixing noise and BBSAN for the conditions examined, but 
some discrepancies remain that are due to the accuracy of the steady RANS turbulence 
model closure, the equivalent sources, and the use of a simplified vector Green’s function 
solver of the linearized Euler equations. 
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Constant associated with the dilatation term 

Coefficients of anisotropic turbulence 

Constant associated with the unsteady forcing term 

Speed of sound 

Coefficients of integral scales 

Nozzle exit diameter 

Fully expanded jet diameter 

Temporal and spatial variation of velocity cross-correlation 
Frequency 

Unsteady force per unit mass 

Unsteady force per unit mass associated with fluctuating shock and turbulent velocities 
Turbulent kinetic energy 

Turbulent kinetic energy associated with fine-scale turbulence 
Wavenumber of the shock pressure in the axial direction 
Integral length scale 

Magnitude of the cross-stream length scale 

Mach number 

Design Mach number 

Fully expanded Mach number 

Pressure 

Pressure due to the shock-cell structure 

Logarithm of the pertubation pressure normalized by 7 

Spectral density 

Strouhal number 

Lighthill stress tensor 

Time 
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Time averaged velocity vector 

Fully expanded jet velocity 

Integral velocity scale 

Spatial coordinate or observer position 

Source location 

Ratio of specific heats 

Dissipation 

Angle from the jet inlet axis to the observer 

Logarithm of the non-dimensional acoustic pressure divided by 7 

Periodic vector Green’s function of the LEE 

Retarded time 

Integral time scale 

Shear stress 

Dilatation source term or angle from the downstream jet axis 

Source term of Lilley’s equation 

Density 

Separation vector, (£,f7,C)> between two source points 
Specific dissipation of turbulent kinetic energy 
Radial frequency 

Broadband shock-associated noise 
Computational fluid dynamics 
Linearized Euler equations 
Nozzle Pressure Ratio 
Power spectral density 
Reynolds- averaged Navier-Stokes 
Sound pressure level 
Shear stress transport 
Total Temperature Ratio 


Introduction 

This paper presents an acoustic analogy for jet noise that is based on the Euler equations. The model 
contains equivalent sources for the turbulent mixing noise and the broadband shock-associated noise (BB- 
SAN). The comprehensive model development is based on two separate acoustic analogies, one for the mixing 
noise and the other for the BBSAN. The evaluation of the comprehensive analogy or the summation of the 
spectral density of the two independent analogies results in a prediction for total jet noise. These models 
are designed for a wide range of Mach numbers and temperature ratios. In the context of this theory there 
are no assumptions of shear-, self-, large-scale, or fine-scale noise sources or spectra. 

The acoustic analogy for turbulent mixing noise is inspired by the work of Morris and Farassat, 1 Morris 
and Boluriaan, 2 Tam and Auriault, 3 and physical measurements. The equivalent sources of the acoustic 
analogy consist of the dilatation and unsteady forces per unit mass. The dilatation term is consistent 
with that of Tam and Auriault 3 and the unsteady forces per unit mass are equivalent to those of Morris 
and Boluriaan. 2 These choices of terms allow for the acoustic intensity to scale with increasing jet Mach 
number and temperature. The acoustic analogy for BBSAN is the model of Morris and Miller. 4 A statistical 
model of the two-point cross-correlation of the turbulent velocity fluctuations describes the properties of the 
turbulence. Frequency dependence of the integral length scale is based on the work of Lieb and Goldstein 5 
and Morris and Boluriaan. 2 

The models are dependent on the vector Green’s function of the linearized Euler equations (LEE) and 
a steady Reynolds-Averaged Navier-Stokes (RANS) solution. Steady RANS represents a compressible base 
flow of the LEE. The computational domain contains the complete nozzle geometry and fully resolves the 
jet plume. Steady RANS solutions are performed with the NASA Langley Fully-Unstructured Navier-Stokes 
(FUN3D) Computational Fluid Dynamics (CFD) package. The numerical solutions depend on the closure by 
the Menter 6 Shear Stress Transport (SST) turbulence model, the nozzle geometry, and boundary conditions. 
Additional details regarding the solution methodology and computational domains are discussed in a later 
section. Predictions are then compared with measured results. 

Characteristics of Jet Noise 

Figure 1 shows normalized power spectral densities (PSD) of pressure from a supersonic off-design jet in the 
far-held, in which the dominant noise types are labeled. The x-axis represents non-dimensional frequency in 
terms of Strouhal number, St = fDj/uj , where / is the frequency, Uj is the fully expanded jet velocity, and 
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Dj is the fully expanded diameter. The y-axis represents lossless Sound Pressure Level (SPL) per unit St, 
which is PSD plus 101og 10 [uj/Dj]. ip is the angle from the jet upstream axis to the observer. Mixing noise, 
labeled in Fig. 1, is induced by the turbulent mixing in the jet, and is dominant in the downstream direction 
of the jet flow. BBSAN is labeled in Fig. 1 and is dominant at the sideline and upstream direction at mid 
to high frequencies. Discrete components can sometimes be observed in the spectra if the jet is operating 
off-design and are called screech. 

Mixing noise occurs at all jet Mach numbers and temperature regimes. At sufficiently high jet Mach 
numbers the directivity resembles a cardioid. The sound directivity is due to a combination of the nature 
of the source and the sound propagation through the jet itself. Tam 7 showed how turbulence contributes to 
far-held jet noise by differentiating the source with hue- and large-scale noise similarity spectra. Jet noise 
has been reviewed by many investigators, including Ffowcs Williams 8,9 and Goldstein, 10 and specifically for 
supersonic jet noise, by Tam. 11 When the convective Mach number of the jet is high enough another mixing 
noise mechanism is enabled that is called Mach wave radiation. For an overview of Mach wave radiation see 
Tam. 12 

A pressure mismatch between the jet how and the ambient environment creates a shock cell structure as 
shown in the schlieren of Fig 2. Shock-associated noise is divided into two parts called screech and BBSAN. 
Screech, hrst observed and described by Powell, 13 consists of discrete tones and is not always observed in 
the far-held PSD. It is caused by a feedback loop consisting of acoustic waves propagating upstream due to 
their reflection by shock waves and consequently being rchected downstream by the airframe or nozzle. For 
more information on screech see the comprehensive review by Raman. 14 BBSAN consists of multiple broad 
lobes and is dominant over a wide range of frequencies in the sideline and upstream directions relative to the 
jet how. It occurs due to large-scale coherent turbulence interacting with the shock waves in the shear layer. 
The resultant acoustic waves combine constructively and destructively at observer locations to create the 
broad lobes. It was hrst characterized experimentally and modeled by Harper-Bourne and Fisher. 15 Notable 
models for the prediction of BBSAN are by Harper-Bourne and Fisher, 15 Tam 16 Tam, 17 and Morris and 
Miller. 4 


Previous Approaches and their Relation to the Present 

Lighthill, in his now famous works on general theory 18 and turbulence as a source of sound, 19 introduced 
the acoustic analogy. By combining the time derivative of the continuity equation with the divergence of the 
momentum equation, an inhomogeneous wave equation results, 

d2 P 2 9 2 p _ 

dt 2 °° dxidxi dxidxj 

where T,y = puruj — Tij + (p — pc^jdij is the Lighthill stress tensor. A solution for p can immediately be 
formed by the convolution integral of the Green’s function of the wave equation with the double divergence 
of Tij. Lighthill constructed a model for Tij consisting of convecting quadruples that resulted in an estimate 
of jet noise intensity scaling as M 8 . 

Proudman, 20 created an expression for the acoustic power from isotropic turbulence using both dimen- 
sional analysis and the acoustic analogy approach of Lighthill. 18 Proudman showed that the acoustic power 
per unit mass from isotropic turbulence is proportional to aeM 5 , where e is the dissipation, M is the Mach 
number based on turbulent velocity fluctuations, and a is a constant determined from the statistical theory 
of Heisenberg. 21 

Curie 22 examined the effects of noise from surfaces by extending Lighthill’s 18 theory. Later, Ffowcs 
Williams and Hawkings 23 constructed generalized forms of the continuity and momentum equations in- 
volving two regions of the flow-field and created an acoustic analogy that involved the generalized density 
perturbation, a generalized Lighthill stress tensor, and two additional source terms that governed the sound 
from turbulence and the effect of moving surfaces. The latter two terms have characteristics associated with 
traditional acoustic dipole and monopole theory. 

Perhaps the most important development in the field of aeroacoustics after the work of Lighthill 18,19 was 
by Lilley 24 who created what is now known as Lilley’s equation. Lilley’s equation is, 


£>V _ D , 2 2 ,, _ dc 2 D dr ' 2 dui <9V 

Dt 3 Dt, ^ dy -2 Dt dy 2 + dy 2 dy\dy 2 


A (y,t) 


(2) 
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where r' = (l/ 7 )ln p, D / Dt = d/dt + u\d/dx \ , and Ui and c are the mean velocity and speed of sound 
respectively. A is a nonlinear source term that involves the base flow and perturbation quantities and 
doesn’t contain r' . This represents a major advantage over Lighthill’s approach. Also, the operator contains 
gradient terms involving the averaged velocity and mean enthalpy. This point is stressed as the effects of 
sound refraction and propagation due to the jet shear layer are taken into account directly and are separated 
explicitly from the source model. Unlike the approaches of Lighthill, Proudman, or Curie, we follow Lilley’s 
approach in the spirit of not having dependent variables on the right hand side of our equations when the 
acoustic analogy is formed and include refraction effects in the operator. 

Khavaran 25 used a prediction methodology that followed the ‘MGB’ code of Mani et al. 2R but differed 
significantly by informing the equivalent sources with a steady RANS solution. This is unlike the approach 
of Mani et al. 26 which calculated the jet flow following the approach of Reichardt. 27 The approach of 
Khavaran 25 has the benefit of including the nozzle geometry in the calculation and shock waves within the 
jet plume. In this work we follow Khavaran 25 by including the nozzle geometry and shock waves with steady 
RANS solutions. 

Goldstein 28 showed that the Navier-Stokes equations can be rearranged so that the right hand side 
consists of sources that are dependent on generalized external stresses and energy flux perturbations. The 
left hand side operator consists of the linearized Navier-Stokes equations. The strength of the approach can 
be demonstrated by using different base flows. As an example, a parallel flow assumption for the base flow 
results in the inhomogeneous Pridmore-Brown equation (like Lilley’s 24 equation). 

Tam and Auriault 3 developed a mathematical theory to predict the noise from fine-scale turbulence based 
on the premise that jet mixing noise consists of two distinct spectra: fine-scale and large-scale, based on 
experimental observation of mixing noise developed by Tam et al . ' Tam and Auriault 3 used a gas kinetic 
theory to relate the time rate of change of the turbulent kinetic energy of the fine-scale turbulence per unit 
mass in the moving frame to the source strength. They argued that the acoustic pressure is related to the 
turbulent kinetic energy of the fine-scale turbulence as, p = (l/3)p < u 2 >= (2/3 )pK s , where u are the 
fluctuating quantities of velocity and K s is the turbulent kinetic energy of the fine-scale turbulence. The 
governing equations are the LEE with the right hand source of the momentum equation equal to —dq s /dxi. 
A closed form solution for the acoustic pressure is found by the convolution integral of the vector Green’s 
function of the LEE with the source term. The predictions of the fine-scale mixing noise agree very well with 
measurement and scale correctly with both jet velocity and temperature. We elect to use a source term for 
the fluid dilation that is equivalent to Tam and Auriault. 3 

Morris and Farassat 1 compared the acoustic analogy approach of Lighthill, 18 the model for fine-scale 
turbulence of Tam and Auriault, and a modified form of Lighthill’s approach based on the same assumptions 
of Tam and Auriault. It was shown that with consistent descriptions of the statistics of turbulence within 
the jet plume, equivalent prediction formulas at the sideline direction of the jet can be formed based on 
the approach of the acoustic analogy (or Tam and Auriault). The work of Morris and Farassat 1 inspired an 
acoustic analogy to be developed by Morris and Boluriaan 2 that is based on the Euler equations, and consists 
of the acoustic pressure written in terms of its logarithm (following Lilley 24 ) and the acoustic velocity, but 
retains viscosity when determining the turbulent mean flow. A left hand side operator of the LEE is formed 
with equivalent sources of the dilatation and an unsteady force per unit mass. The approach within this 
paper draws heavily upon Morris and Boluriaan by using the same source term for the unsteady forces per 
unit mass and retains the advantage of formulating the operator of the LEE with only four unknowns. 

Harper-Bourne and Fisher 15 proposed that BBSAN depends on the coherent interaction of turbulence 
in the jet shear layer with the nearly periodic shock cell structure. This is modeled as a series of correlated 
point sources that radiate constructively or destructively. The prediction scheme depends on knowledge of 
the rate of decay of the turbulence correlation between shocks, as well as the characteristic spectral shape 
of the radiated noise generated by each interaction. Tam 16 developed a stochastic methodology for the 
prediction of BBSAN. It depends on the large-scale turbulence of the jet shear layer being modeled as a 
superposition of random instability waves interacting with the nearly periodic shock cell structure. This 
interaction yields an interference pattern of traveling waves with phase velocities higher than the speed 
of sound. The instability waves give rise to noise radiation at large angles from the jet downstream axis. 
Tam 1 ' modified the method of Tam 16 to include moderately heated jets and additional corrections for over- 
or under-expanded jets. Neither of the approaches of Harper-Bourne and Fisher 15 or Tam 16 are based on 
the acoustic analogy and they are incompatible mathematically with the previously discussed approaches. 
Morris and Miller 29 developed an acoustic analogy for BBSAN based on the Euler equations, and it drew 
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heavily on the work of Harper-Bourne and Fisher 15 and Tam and Auriault. Their approach is used here as 
it is compatible with the current mathematical formulations with minimal modification. 


Mathematical Theory 

The acoustic analogy presented herein is based on that of Morris and Farassat. 1 The governing equations 
are the Euler equations, 


Dir dut 
Dt ^ dxi 


( 3 ) 


Dui 2 ^ 7r 
Dt C dxi 


= 0 


( 4 ) 


where D/Dt is the material derivative and ir = 'y~ 1 ln(p/p oc ). By linearizing the governing equations about 
a mean flow and retaining terms on the left hand side that are linear in fluctuation, then the inhomogeneous 
LEE are formed. 


dir' _ dir' du ' 

dt J dxj dx-i 

and, 


( 5 ) 


dK 

dt 


_ du 1 , . dui 

"'dXj ' " J d.r j 


,2 dn 

dxi 


fi 


( 6 ) 


where 9 is the dilatation rate and fi is the unsteady force per unit mass. In both the equivalent sources 
only the second order fluctuations are considered to contribute to the sound field. Note that it' ~ + 

p'/Poo) ~ 7 _1 pVPoo- The vector Green’s function of the LEE is defined as, 


and, 


DoK dU 9i 

Dt dxi 


5{x - y)S(t - T)S 0n 


( 7 ) 


D ° u gi . n du i ,-2 dn 9 x( w , u 

-m +u «W i +c W t =s{x - y)sif - T)5> '" 

where D a is the material derivative about the base flow. Let the vector Green’s function be periodic, 


(8) 


1 f°° 

Kg(x,t\y,T) = — J n™(x\y;u})exp[vx(t - r)}du> 

/♦OO 

Kg(x\y,u)= / iTg(x,t\y,T) exp[—iu(t — r)]dr. 


(9) 

(10) 


J — OO 

It is now simple to write the fluctuating far-held pressure as a convolution integral of the vector Green’s 
function and the equivalent sources, 


/ oo ro o 3 

... / TV°(x,t\y,T)B(y,T) + ^ t|j/, r)/„(y, r)drdy. (11) 

-oo J- OO n=1 

The spectral density is dehned as the Fourier transform of the autocorrelation of the fluctuating pressure, 

/ OO 

p’(x, t)p'(x , t + r) exp [— iwr]dr. (12) 

-OO 


Equations 11 and 12 yield, 
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Cj) Poo C o 


(x, y , -cj)7T°(x, y + r/, o;)0(y, r)0(?/ + y,t + r) 


+EE 7Tg (®, y, (®, y + y, uj)fn(y , T)f m (y + y, t + r) > exp[-*wr]drrfydy. 

n=l ra— 1 J 

By definition the periodic vector Green’s function of the LEE is, 

Kg*{x,y,u>) = TTg(x,y, -oj) (14) 

and assuming that the observer is in the far-field relative to the source, we relate two closely placed source 
points in the jet, 

n™(x,y + y,w) ~ 7r™(a;,y,w)exp — ^(xi£ + x 2 y + x 3 () . (15) 

|_ Coo X 

Using Eqns. 14 and 15 with Eqn. 13 yields, 


$ (x, Poo C o 


|7T*°(x, y, u)i T° g (x, y, w)0(y, r)6(y + y, t + r) 


+ EE 7T* n (x, y, u)n™(x, y, w)/ n (y, r)/ m (y + y, t + r) 


x exp (xi£ + x 2 y + X 3 C) exp [— zwr] drdr/dy. 

|_ Coo J 

Following the source model of Morris and Boluriaan 2 (which is equivalent to Tam and Auriault 3 ) for the 
dilatation rate, 

0(y,T)9(y + ri,t + r) = (17) 

' s 

where A s is a constant. The source model of Morris and Boluriaan 2 for the unsteady force per unit mass is, 


fn{y,r)f m (y + rj,t + T) = gz K/c -°°) u s 


where B s is a constant. A model must be formed for E(r), r). 

Following Ribner’s 30 postulate, the spatial and temporal terms of the model are separable, 

. . r-lrll r-(e-Ur) 2 l f — (?7 — 1 [ — (£ — urr) 2 "| , , 

E(n h t) = exp exp — x— exp exp — X— (19) 

I T s \ L J L l y \ L l i J 

The integrations involving y and r can be performed analytically. These integrals and associated variables 
are isolated, 


— |r|l M^-ur) 2 ] r-(y-ur) 2 l [ — (£ — wt ) 2 

exp exp exp exp 

> ^ J L g J L s J L 


x exp (xi£ + X 2 y + X 3 Q exp \—iu>T] d^dyd^dr 

CooX 


The integrations with respect to £, Q, y, and r, results in, 


2'7T 3 ^ 2 C^ 0 l x lylzX 2 T s 


C %>X 2 + (c a C X + UX 1 + UX 2 + WX 3 ) 2 T 2 L 0 2 


(llx\ + l 2 v xl + Z 2 x§) w 2 


4c 2 x 2 


The integrations of y and r require the assumptions r s > 0, l t > 0, > 0, x > 0, where each of these 

is real. These are reasonable assumptions physically as the turbulent kinetic energy at the boundary is zero 
and we restrict our modeling of the jet flow to x > 0, where x = 0 and —D/2 < r < D/2 represents the 
nozzle exit. We can now write the spectral density as, 
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/ oo /*oo /»oo 

/ / 

-oo J — oo J — 


2ir 3/2 cl 0 l x lyl z x 2 T s 


/_ oo + (CooO: + UX1 + + wx 3 ) 2 t 2 lo 2 

/ / \ 4 3 3 

l2 ( w «/Coo) 


X ^ n g 0 ( x i Vi u)-kI{x, y , iu)A 2 V “ s/ J 2 °° ; + H 


2 (u s /c 00 ) 2 u 2 


ll 


( 22 ) 


x exp 


(^i + + ^s§) a; 2 

4c^,a: 2 


dy. 


Broadband Shock- Associated Noise 

A separate approach has been developed to predict the BBS AN component of the total jet noise. This 
approach is based on an acoustic analogy that follows a similar methodology to that outlined for the mixing 
noise. Thus, the two methodologies are compatible and the results can be summed. The BBSAN component 
of the jet noise is calculated using the acoustic analogy of Morris and Miller 4 and its development is briefly 
discussed here. 

The equivalent source terms depend on the product of the jet’s turbulent velocity fluctuations and the 
shock cell structure pressure, and a statistical model is used to describe the properties of the jet turbulence. 
The model of Morris and Miller is, 


S ( x , to) = 


IQn^c^x 2 J^ c 


Kl 2 

~r L Ps(y)Ps(ki,y2,y3 ) 

IT a 


uj 2 t 2 exp —l 2 (k i — wcos0/coo) 2 /4 — u 2 l\ sin 2 0/4 c^, 

1 + (1 — M c cos 8 + u c ki/uj) 2 u> 2 T 2 dkidy 


(23) 


where M c is the local convective Mach number of the turbulence, Ps{k\,y 2 ,y 3 ) is the Fourier transform of 
the shock pressure (p — p^) in the axial direction, k\ is the wavenumber component in the axial direction, 
and u c is the convective velocity of the turbulence in the axial direction. 

This BBSAN model uses the vector Green’s function of the LEE in a quiescent environment for propaga- 
tion. Realistic propagation effects can be taken into account by using the approach described by Miller and 
Morris 31 or Henry et al 32 Since the propagation of the BBSAN through the shear layer is not being captured 
by the choice of the Green’s function, an empirical correction suggested by Tam 17 was used to calculate the 
scaling of BBSAN with increasing temperature, but could more accurately be taken into account as shown 
by Miller. 33 

The integral scales of turbulence in Eqns. 22 and 23 are both based on steady RANS solutions. The 
length scales within the BBSAN model are chosen to be independent of frequency as there is little effect 
on the predicted BBSAN with frequency dependent length scales. Since BBSAN is caused by large-scale 
turbulent structures, small-scale turbulence is not relevant to the calculation. 


A Partially Comprehensive Model 

The models for spectral density of mixing noise and BBSAN are now combined into a single model equation. 

/ OO nOO 

••• / {?r *g O (x,y,oj)n 0 g (x,y,u>)8(y,T)8(y + ri,t + T ) 

-oo J — OO 
3 3 

+ Y^ 7r s"( a; ’y’ a; ) 7r r( a: ’y’ w )/™(y’ r )/ m (y + 7 ?’ t + r ) 

n— 1 m = 1 
3 3 

+ Y1 Y, 7r *g n ( x ^y’ UJ ) 7r T( x ’y^)fn(y^ r )fm(y + V: t + T )} 


x exp 


-(zi£ + x 2 y + x 3 () 


exp [-iwr] drdridy. 


where, 


( 24 ) 
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and u' s j and v! ti 
and Miller, 29 


du’u 


du' 


= <*> 

are the velocity perturbation due to the shocks and turbulence respectively. Following Morris 


fn(y^ T )fm(y + r i,t + r) = 


timn (■ y)Ps{y)Ps{y + ri)k{y)E(ri, r) 


p{y) 2 u 2 (y)iUy) 

Substituting in the source terms from the mixing and BBSAN model, 


(26) 


/ oo poo r 

... I Tr* g 0 (x,y,u>)Tr° g (x,y,u)A 2 s 
-oo J — OO L 


(u 8 / ^oo) 


3 3 


, *n/ \ m( \ I d 2 ( W s/ C oo) u s , a mnk{y)ps{y) I ~ n \ r * 7 1 

+ 2^ 1^*9 ( x ^y^) n g [ B S J 2 + 0 ^-2-2y 2 — / Ps{ki,y 2 ,ys)exp[-ik 1 y 1 \ 


xE (77, r) exp 


2?r p 2 u 2 ll j_ oo 

—iuj 


CocX 


-(x 1 £ + x 2 'n + x 3 () 


exp [— icor] dTdr/dy. 

(27) 


where, 


/ OO 

Ps{y) exp [ik\y\\dyi 

-OO 


and the prefactor (27 t) 1 is retained in the forward transform. Integrals of Eqn. 27 are isolated, 


POO poo 


E ( rj , r) exp 


' — oo */ — OO 


-{xi£ + x 2 ri + x 3 () 


exp [— iur] drdr) 


(28) 


(29) 


The evaluation of expression 29 and the combination of Eqn. 27, after simplification, results in, 


/ OO pOO p 

/ / 

-OO J — OO J — OO '"'OO 


2’K 3/2 cl 0 l x lyl z X 2 T s 


C^X 2 + (■ UXi + VX 2 + wx 3 + CooXYt 2 ^ 2 


I *0/ \ 0 / \ ,2 ( M s/ c oo) 

X 1 7T g u (x,y,u>)7r“(x,y,uj)A s ^ 


+ EE 7r g n ( x ’y’ u ') 7r a l ( x >y’ u ') ( B ‘> 

n—1 m = 1 

C ’ mn k(y)p s (y) 


( u a /coo) 2 u 2 s 

ll 


(30) 


27 Tp 2 U 2 l^ 


Ps(ki,V 2 , 2/3) exp[— ifcij/i] 


x exp 


-( l l X l + l v X 2 + l z X t)° 


V 

4c 2 x 2 


dy. 


Equation 30 is a partially comprehensive analogy that contains noise sources for all mixing noise compo- 
nents of jet noise and the broadband component of the shock-associated noise. It does not contain a source 
term for the discrete component of shock-associated noise. It also assumes there is no additional noise source 
due to the interaction of the mixing noise source with the shock-associated noise source. This would require 
a model for an equivalent source of the form f n (y, t)f^(y + r/, t + r) or perhaps 9(y, r)f^(y + rj,t + t). At 
this time it is unknown what models are appropriate for these potential source terms. 


Frequency Dependent Integral Scales 

Equations 22, 23, and 30 require the integral scales of turbulence. Jet flows exhibit self similarity and 
the length scale grows linearly with increasing streamwise distance. An empirical model of the growth of 
the length scale with streamwise distance might be approximated as, l x (yi) — 0.138j/iD, where y\ is the 
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streamwise distance from the nozzle exit to the source. Here, the integral scale of turbulence is found from 
K and e directly from the steady RANS solution as, 

lx{v) = cimi x K(y)/e(y) (31) 

which is the simplest approach that retains the correct dimensions. 

Morris and Zaman 34 examined the variation of the length scale with St. They observed that from low to 
mid frequencies the integral length scale is relatively independent of St and at high frequencies the integral 
length scale falls off as the inverse of St. Morris and Boluriaan 2 incorporated frequency dependent length 
scales in their model for jet noise and demonstrated improved predictions. Lieb and Goldstein 5 examined 
frequency dependent length scales and showed their effect on jet noise predictions. After examining multiple 
choices for an appropriate frequency dependent length scale we select, 

i x (y, st) = Uy) l ~ e ^st fSt] ( 32 ) 

where c/ = 11.25 is a coefficient that is set to match the observations of Morris and Zaman. 34 The temporal 
and velocity scales must also be related to the steady RANS solution. A simple dimensional model that 
follows Eqn. 31 is selected. The temporal and velocity integral scales of turbulence are, 

r s (y) = Cr mix K(y)/e(y) (33) 

and, 

u s {y) = u amix y/2K(y)/3 (34) 

where c Tm i x and u am i x are constant coefficients. 

The constants shown in Eqns. 30, 31, 33, and 34 are specified by performing a parametric study relative 
to a single jet condition and observer angle with various sets of experimental data. The coefficients are 
A s = 0.25, B s = 2.51, c T mi X = 0.06, u am i X = 1.00, and c; TO j x = 0.40. The cross-stream length scales are 1/3 
of the streamwise length scale, which is a typically accepted value. The effect of using a frequency dependent 
length scale for mixing noise is justified in terms of representing the physics of turbulence in jets and has 
improved predictions as shown by Lieb and Goldstein. 5 

Steady RANS Solutions 

The steady RANS solutions are found using the FUN3D solver developed at NASA Langley Research 
Center. For more information on FUN3D see Anderson and Bonhaus 35 or Nielsen. 36 The BBSAN model is 
compatible with other steady RANS solutions and was demonstrated with the Wind-US code (see Nelson 37 
for details on Wind-US). 

In this paper multiple axisymmetric nozzles are examined. A single plane in the streamwise and radial 
direction of the nozzle is constructed that contains the nozzle interior, exterior, and plume region. This 
domain contains the exact axisymmetric geometry used in the corresponding experiment. The axisymmetric 
plane is rotated 90 deg. in the azimuthal direction about the centerline axis of the jet (a: = 0 and y = 0). 
25 planes are created in the azimuthal direction. This operation yields a three-dimensional grid consisting 
of hexahedrons except on the centerline axis that consists of prisms. 

Singlestream axisymmetric nozzles of this study are the SMC000, SMC014, SMC016, and SMC018 of 
NASA Glenn Research Center. They all have exit diameters of 0.0508 m. Figures 3 and 4 show a region 
of the computational domain of the SMC000 and SMC016 nozzles respectively. The domain extends 1001? 
downstream from the nozzle exit, 501? in the radial direction from the centerline, and 5 D upstream from 
the nozzle inlet. Another nozzle examined is the Boeing Company axisymmetric dual-stream nozzle shown 
in Fig. 5, that has a primary exit diameter of 0.06223 m and a secondary exit diameter of 0.152 m. The 
ranges of the dual-stream nozzle domain are the same as the single-stream domains. Each of these figures 
are normalized about the exit diameter of the core stream. 

Numerical solutions are governed by the steady RANS equations closed by the Menter SST 6 turbulence 
model. Roe 38 flux vector construction is used for spatial discretization and is second order accurate. A Roe 39 
‘modmin’ flux limiter is used in cases that fail to converge due to very large fluxes in the initial transient 
solution. This approach results in a total variational diminishing scheme. An ‘inlet 1 boundary condition 
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is set at the nozzle inflow that holds total pressure and total temperature constant. At the nozzle wall an 
adiabatic no-slip condition is set and the normal flux is zero. Far-field boundaries are set as ‘freestream’ with 
A/oo = for stability and convergence. An ‘outflow’ boundary condition is set at the downstream faces 
that hold the static pressure equal to the ambient, p^ = 101325 Pa. Finally, the x — y and x — z planes, 
that intersect the nozzle centerline axis, are set as symmetric boundaries. These boundary conditions are 
the same as those used by Miller and Veltin. 40 

The flowfield is initialized as quiescent at the start of the simulation. Initial iterations use a Courant- 
Friedrichs-Lewy (CFL) number of 1.0 and increase to a CFL number of 50.0 after 250 iterations. The 
total number of iterations required for convergence varies depending on the nozzle and boundary conditions. 
Typically the total number of iterations required for convergence is approximately 8000. Figure 6 shows the 
global residual per iteration for a typical case. The global residual of p has been reduced by three orders of 
magnitude and is the smallest overall reduction. In all simulations the domain is decomposed using parallel 
graph partitioning and fill-reducing matrix ordering and run on multiple processors. When the solver finishes 
the solution is inspected visually, and the residuals are examined to determine if the case has converged. The 
solutions consisting of p, u, v, w, p, K, and f l are written to file. The validation of the FUN3D steady RANS 
solutions for single-stream jets was performed and is presented in an accompanying paper of this conference. 

Figure 7 shows example steady RANS solutions in the vicinity of the nozzle exit for two cases. The 
conditions are a) Afo = 1.00, Mj = 1.47, and TTR = 3.20 (Total Temperature Ratio) and b) M r j ; = 1.50, 
Mj = 1.50, and TTR = 1.00. Single planes are shown with contours of the streamwise velocity component. 
Clearly, the nozzle geometries are captured within the simulation. The first case shows that the shock cell 
structure is captured if the jet is operating off-design and is taken into account both in the steady RANS 
simulations and the predictions for mixing noise and BBSAN. 

Implementation 

Equations 22, 23, and 30 are implemented in a computer program. The program uses structured or 
unstructured steady RANS solutions and converts the specific dissipation, f l, to dissipation e = 0.09A'f2. 
The steady RANS solution is interpolated using the inverse weighted distance method onto a new structured 
grid that encompasses the jet plume. Numerical integration is performed using the new structured domain 
to approximate the integrals of the model equations. Evaluation of the shock-associated source term requires 
p s . Additional details for calculating p s are shown in Morris and Miller. 29 The source term for BBSAN has a 
coefficient a mn , and term 7r™, that are approximated jointly with the assumption of Proudman. 20 Integration 
ranges and grid point density are varied until the solution is grid independent. 

Tv" is governed by Eqns. 7 and 8. Unfortunately, a full numerical solver for n™ is not available to 
the author but numerical solution techniques have been demonstrated by Tam and Auriault. 41 Here, a 
simplified approach has been undertaken. The spreading rates of jets are relatively small and a locally 
parallel flow assumption is made. Equations 7 and 8 are written about a locally parallel base-flow at each 
streamwise location. They are reduced to a single third-order differential equation cast in an adjoint form. 
The steady RANS solution is used to provide the spatial derivatives of density, streamwise velocity, and 
pressure. Numerical solutions of the differential equation yield the Green’s function of Lilley’s equation and 
are related to tt™. Full details of the approach to calculate 7r™ are described by Miller and Morris. 42 Using 
a locally parallel mean flow will have some effect on the predictions in the downstream direction of the jet, 
particularly for noise sources located within the turning point of the flow. These have been quantified by 
Tam and Auriault. 41 The parallel flow assumption does not affect the model development presented in the 
previous section. 


Results 

The developed model is assessed over a wide range of jet Mach numbers and temperature ratios. Single- 
stream jet cases are shown in Table 1. Each row represents a case with associated design Mach number, fully 
expanded Mach number, and total temperature ratio. The exit diameter of the single-stream cases (SHJAR 
database) is 0.0508 m. Dual-stream jet cases are shown in Table 2. A subscript p represents a condition 
associated with the core stream and a subscript s represents a condition associated with the bypass stream. 

The first set of predictions are shown in Figs. 8 through 10 and represent jet flows from a convergent nozzle 
with Mj > 1. The nozzle geometry is the SMC000 convergent nozzle of Bridges and Brown. 43 The observers 
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are located at six angles from the nozzle exit at a distance of R/D = 100. Measured data from Bridges and 
Brown are shown in red. Predictions of mixing, BBSAN, and total noise are shown as blue, green, and black 
lines respectively. Generally, in the upstream and sideline direction the predicted BBSAN dominates the 
mixing noise from mid to high frequencies and is dominated by mixing noise in the downstream direction. 
This is consistent with experimental observation and the breakdown of jet noise spectra by Viswanathan. 44 

Figure 8 shows spectra from a Mj = 1.47 and TTR = 1.00 jet. The measured data clearly shows 
BBSAN and screech. At the sideline and upstream directions (ip = 50 through 90 deg.) the mixing noise 
under-predicts the low frequencies in the range of St = 0.02 through 0.2. In the rest of the frequency range 
BBSAN is dominant. The peak BBSAN frequency is fairly accurately predicted by the model at all angles, 
and the amplitude is in particularly good agreement, even with strong screech tones in the sideline and 
upstream directions. At higher frequencies (St = 8 and above) the BBSAN model often under-predicts the 
measurement by 2 to 4 dB. In the peak noise direction, ip = 130 deg., or near the jet axis, ip = 150 deg., the 
mixing noise is about equal to or dominates the BBSAN respectively. 

A second case is shown in Fig. 9 with Mj = 1.00 and TTR = 2.70, that represents a sonic heated 
jet. Since the jet is sonic it can not support shock waves and the BBSAN source in the model equation is 
zero. The measured spectra do not show any shock-associated noise component. The predicted mixing noise 
magnitudes agree within 4 dB of the measurement at all angles except ip = 130 deg. where the prediction 
is approximately 6 dB below measurement. At the sideline direction the mixing noise prediction shape is 
similar to measurement but over-predicts the peak frequency by St ~ 0.25. This is not the case in the 
upstream direction at ip = 50 deg. where the peak frequencies agree. Near the jet downstream axis the 
effect of the ‘large-scale’ source is well captured in magnitude. At very high frequencies in this direction the 
model under-predicts the measurement by 4 dB near St = 10. 

Figure 10 shows the spectra for a Mj = 1.24 and TTR = 3.20 jet. This flow is weakly off-design relative 
to that shown in Fig. 8 (observe the relative strengths of the BBSAN) and has additional heating. The 
relative magnitudes of the mixing noise compared to the measurement show similar trends to the previous 
case even with the addition of heating and a slight increase in Mj. The frequency peaks of BBSAN are 
aligned with measurement but under-predict the magnitude, especially in the downstream direction. At 
some locations such as the sideline, it is illustrative to observe that the total noise is a combination of both 
mixing and BBSAN, especially so near St ~ 1.0. 

The predictions shown in Figs. 8 through 10 are representative of the first nine cases shown in Table 1. 
Overall, the amplification of mixing noise due to the propagation through the jet shear layer is captured. 
This is the case even with the locally parallel mean flow assumption. 

We now examine subsonic jet predictions from the convergent SMC000 nozzle. These are shown in 
Figs. 11 through 13. In all these cases the source of BBSAN is zero. Figure 11 shows a high speed heated 
subsonic jet operating at Mj = 0.90 and TTR = 3.20. The predicted peak magnitudes are within 2 dB of 
the measurement except for ip = 130 deg. where the prediction is 7 dB lower than measurement, but the 
evolution of the spectral shapes from the upstream through downstream direction show a positive trend with 
respect to measurement. As an example, the predictions at ip = 50 relative to ip = 150 deg., have a much 
sharper high magnitude peak. There is also noticeable error in the prediction of the peak frequency that 
varies between —0.10 and 0.3 St. 

Figure 12 shows a prediction for a Mj = 0.70 and TTR = 1.00 jet. This case, which is an unheated flow, 
shows similar trends relative to the heated Mj = 0.90 case, implying that the temperature scaling of the 
model is reasonable at subsonic speeds. Some differences can be seen relative to Fig. 11. For example, at 
ip = 150 deg. the lower speed case under-predicts the peak jet noise by 4 dB but predicts the peak magnitude 
accurately at ip = 130 deg. Unfortunately, the peak frequencies show some of the same shifting as seen in 
the supersonic cases. The absolute magnitude of the predictions relative to measurement at the peak of the 
spectra differ by ± 3.5 dB. 

A final subsonic prediction is shown in Fig. 13 for a Mj = 0.50 and TTR = 1.00 jet. Similar trends are 
seen relative to the other subsonic cases. The varying Mach numbers and temperature ratios of the subsonic 
jets shown are illustrative of cases ‘SMC00010’ through ‘SMC00017’ of Table 1. 

A series of jets from single-stream convergent-divergent nozzles are now examined. Figure 14 shows the 
spectra from a M^ = 1.185 (SMC014 nozzle) on-design jet operating at TTR = 1.80. This jet, relative to 
the sonic cases discussed previously or the off-design cases, is only slightly more energetic. The spectra are 
nearly on-design, and only extremely weak shocks are observed in the upstream and sideline direction. In 
these directions the prediction agrees very well with measurement from low to mid frequencies but under- 
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predicts the noise at higher frequencies. This is partly due to the shock noise component that dominates 
the noise from St = 0.8 through 10 at the upstream observer angles. At ip = 130 and ip = 150 degs. the 
predicted peak frequency agrees with measurement but the fall-off of the prediction is too large relative to 
measurement . 

We now turn our attention to a family of jets based on the M,j = 1.50 SMC016 nozzle. The predictions 
relative to measurement are shown in Figs. 15 through 18 representing an unheated jet flow ( [TTR = 1.00) 
from the subsonic through under-expanded Mach numbers. Figure 15 represents Mj = 0.784 and the spectra 
should be shockless. Unfortunately, some of the measured data appear to contain very weak shocks even 
though the plume is shock free. 

An over-expanded Mj = 1.294 jet is shown in Fig. 16 and strong BBSAN and screech tones appear. 
The mixing noise predictions at the upstream and sideline directions agree very well with experiment and 
the BBSAN predicts the correct peak frequencies. The BBSAN prediction at the sideline and upstream 
directions, unfortunately, over-predicts the noise near St = 0.4; this causes the total noise prediction to 
be higher than measurement. The prediction at ip = 130 deg. correctly identifies the frequency where 
the mixing noise and BBSAN have equal intensity near St = 1.2. At the highest observer angle the total 
prediction over-predicts the peak frequency by St = 0.15. The BBSAN at this angle has captured the extra 
energy above St = 1 but unfortunately is dominated by an incorrect mixing noise prediction. 

Figure 17 shows the on-design Mj = 1.50 unheated case from the SMC016 nozzle. The measurement 
reflects proper nozzle design as only extremely weak BBSAN can be seen in the upstream direction. Relative 
to the previous case, the predicted peak magnitude and frequencies in the dominant jet noise direction more 
closely match experiments. Unfortunately, at ip = 130 deg. the high frequencies are under-predicted by up 
to 10 dB. 

A final prediction for the M d = 1.50 SMC016 nozzle is shown with Mj = 1.734. In this case the 
BBSAN dominates the spectrum in the upstream and sideline direction and partly in the downstream. The 
measurement contains strong screech tones and these have been shown (Andre et al , 45 ) to have strong effects 
on BBSAN. Except in the upstream direction, the BBSAN peak is not correctly predicted and is often higher 
by St = 0.2 to 0.3. The mixing noise prediction at low frequencies in the sideline and upstream direction 
under-predicts the measured mixing noise. It is likely that for this jet there are at least some non-negligible 
source terms that affect the mixing noise source due to shocks in the flow. 

The last single-stream prediction is performed using the M ( j = 1.80 SMC018 nozzle shown in Fig. 19. The 
jet operates at Mj = 1.80, which is the highest jet Mach number presented (the highest evaluated is 1.828), 
and is moderately heated at TTR = 1.648. The measured spectra are almost shock free. At ip = 130 deg. 
the peak prediction is 2 dB lower than experiment, and at ip = 150 deg. the peak prediction is lower than 
measurement by 4 dB. In the sideline direction the low frequencies are under-predicted. In the dominant 
noise direction the high frequencies are under-predicted. 

A final prediction is shown in Fig. 20 for a dual-stream axisymmetric jet. The core stream operates 
at M d = 1.281 and TTR p = 2.87 and the fan stream operates at Mj S = 0.848 and TTR S = 1.00. The 
nozzle exits of the core and fan streams are off-set as shown in Fig. 5, and both core and fan nozzles are 
convergent. The core diameter, D p , is 0.06223 m and the observers are at R/D p = 100. The core stream 
is operating off-design and the fan stream is subsonic. Since the nozzle exits are off-set, multiple non- 
overlapping integration regions are used to evaluate the model. The measured data in Fig. 20 are courtesy of 
Viswanathan. 46 Both the mixing and BBSAN predictions at high frequencies over-predict the measurement 
by significant amounts at all angles. At lower frequencies (St ~ 0.05) the opposite is true and the predictions 
are lower than measurement, however the measurement does not reflect typical low-frequency trends of jet 
noise. In the upstream and sideline directions the mixing noise and BBSAN peaks are within 2 dB of the 
peak measurement. The quality of predictions of the second case in Table 2 is similar to the first. 

In almost all the predictions shown the noise is lower than measurement at ip = 130 deg. Also, almost 
all predictions at ip = 150 deg. had peak intensity close in magnitude to that of measurement. The 
use of a vector Green’s function solver with a locally parallel mean flow assumption, has shown in these 
calculations that predicting the peak noise intensities is possible, even at very high Mach numbers. The 
directivity and frequencies of the peak intensities in the downstream direction are not always in agreement 
with measurement. With the use of a fully three-dimensional vector Green’s function solver of the LEE it 
could be possible to correct these errors at all Mach numbers. 

The equivalent sources in this work are based on well known models that are related to steady RANS 
solutions. Unfortunately, the accuracy of the steady RANS modeling relative to scaling of the sources with 
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increasing M 3 and TTR is questionable. This is apparent in the assessment paper by Georgiadis et al. 47 
where various models in the same CFD code make vastly different predictions. An assessment of FUN3D 
for the jets examined in this paper had trouble capturing all the flow-features that are essential for accurate 
jet noise prediction for a wide range of conditions. It is hopeful that in the future more accurate steady 
RANS turbulence models will be created. Equivalent sources have substantial potential for improvement. 
One possible improvement is anisotropic effects, such as those examined by Khavaran 48 or Hunter, 49 that 
could substantially improve predictions. 

Spectra varying with increasing observer angle, tp, transition from being dominated by noise from turbu- 
lence that is more universal to that which is dominated by large-scale coherent structures. This transition 
represents, in the theory of Tam et al. 7 the transition between the fine- and large-scale similarity spectra. 
The first two terms of Eqn. 24, in the context of similarity spectra, could be thought of as fine- and large-scale 
sources respectively. The dilatation, unsteady forces per unit mass, and broadband shock-associated sources 
share the same two-point cross-correlation of the turbulent velocity fluctuation. In this theory, there is no 
assumption that the noise source consists of fine- or large-scale similarity spectra and associated source, or 
even traditional shear- and self-noise terms. 


Conclusion 

The Euler equations are rearranged into an acoustic analogy with arguments consisting of the vector 
Green’s function of the LEE, two-point second order cross-correlation of the equivalent sources, and a steady 
RANS solution. The equivalent sources are modeled for jet mixing noise following the work of Tam and 
Auriault 3 and Morris and Boluriaan 2 and for BBSAN following the work of Morris and Miller. 29 A statistical 
model of the two-point cross-correlation of the velocity fluctuations is used to describe the turbulence within 
the jet plume. This statistical source model is shared between the mixing and BBSAN sources. These 
sources reside in a single equation for the total jet noise but do not account for the discrete component of 
shock-associated noise. Propagation effects controlled by the vector Green’s function of the LEE are found 
with a numerical methodology in conjunction with a locally parallel flow assumption. Propagation effects for 
BBSAN are also controlled by the Green’s function of the LEE and are approximated by spherical spreading. 

A single model equation is evaluated for the total noise, mixing noise, and BBSAN with Mj ranging 
from 0.50 to 1.828 and TTR ranging from 1.00 to 3.20. The predictions are compared with experiments of 
jets operating subsonically through supersonically and at unheated and heated temperatures. This includes 
transonic and off-design cases. Predictions generally capture the scaling of both mixing noise and BBSAN 
for the conditions examined, but some discrepancies remain that are due to the accuracy of steady RANS 
turbulence model closures, the equivalent sources, and the use of a simplified vector Green’s function solver of 
the linearized Euler equations. The model makes no assumptions regarding the fine- or large-scale turbulence, 
self- or shear- noise, or convective amplification. Improved predictions can likely be made with the use of a 
fully three-dimensional vector Green’s function solver and improved source and turbulence models. 
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Table 1. Single-stream jet operating conditions. Case names with a **’ are shown in this paper. 


Case 

M d 

Mj 

TTR 

SMC0001 

1.00 

1.00 

1.00 

SMC0002 

1.00 

1.24 

1.00 

SMC0003* 

1.00 

1.47 

1.00 

SMC0004* 

1.00 

1.00 

2.70 

SMC0005 

1.00 

1.24 

2.70 

SMC0006 

1.00 

1.47 

2.70 

SMC0007 

1.00 

1.00 

3.20 

SMC0008* 

1.00 

1.24 

3.20 

SMC0009 

1.00 

1.47 

3.20 

SMC00010 

1.00 

0.70 

3.20 

SMC00011* 

1.00 

0.90 

3.20 

SMC00012* 

1.00 

0.70 

1.00 

SMC00013 

1.00 

0.90 

1.00 

SMC00014* 

1.00 

0.50 

1.00 

SMC00015 

1.00 

0.50 

1.80 

SMC00016 

1.00 

0.50 

2.20 

SMC00017 

1.00 

0.50 

2.70 

SMC0141 

1.185 

1.185 

1.00 

SMC0142* 

1.185 

1.185 

1.80 

SMC0143 

1.185 

1.185 

3.00 

SMC0161* 

1.50 

0.784 

1.00 

SMC0162 

1.50 

1.046 

1.00 

SMC0163* 

1.50 

1.294 

1.00 

SMC0164* 

1.50 

1.500 

1.00 

SMC0165* 

1.50 

1.734 

1.00 

SMC0166 

1.50 

1.828 

1.00 

SMC0181 

1.80 

1.80 

1.00 

SMC0182* 

1.80 

1.80 

1.648 


Table 2. Dual-stream jet operating conditions. Case names with a 4 *’ are shown in this paper. 


Case 

D p (m) 

D s (m) 

M dp 

Mjp 

Mj S 

TTR p 

TTR S 

axi5* 

0.06223 

0.152 

1.00 

1.281 

0.848 

2.87 

1.00 

axi7 

0.06223 

0.152 

1.00 

1.240 

1.358 

2.87 

1.00 
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Figure 1. Sound pressure level per unit Strouhal number resulting from a Md = 1.00, Mj = 1.50, and TTR = 1.00 
jet. The red lines represent the noise that is dominated by turbulent mixing, black lines represent noise that 
is dominated by BBS AN, and green lines represent noise that is dominated by discrete tones (‘screech’). *0 is 
the angle from the upstream jet axis to the observer centered about the nozzle exit at R/D = 100. 
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Figure 2. A schlieren of an over-expanded Mj = 1.60 jet from NASA Glenn Research Center courtesy of 
Zaman et a/. 50 The oblique shock waves and expansion waves are apparent within the jet plume. 



Figure 3. The SMC000 nozzle geometry and a region of the associated computational domain normalized by 
nozzle diameter. 
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Figure 4. The SMC016 nozzle geometry and a region of the associated computational domain normalized by 
nozzle diameter. 


9 




Figure 5. The Boeing dual-stream nozzle geometry and a region of the associated computational domain 
normalized by nozzle primary diameter. 
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Figure 6. Global residuals calculated by FUN3D of the field variables per iteration for the M d = 1.00, M, = 1.47 
TTR = 3.20 jet flow. ’ 




Figure 7 . Example steady RANS solutions of the streamwise velocity component of a) M d = 1.00, Mj = 1.47 
and TTR = 3.20 and b) M d = 1.50, Mj = 1.50, and TTR = 1.00. 
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Figure 8. Comparison of predictions (Case SMC0003) with measurement of the SMCOOO convergent nozzle 
operating at Mj = 1.47 and TTR = 1.00 at a distance of R/D = 100 and angles 
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Figure 9. Comparison of predictions (Case SMC0004) with measurement of the SMCOOO convergent nozzle 
operating at Mj = 1.00 and TTR = 2.70 at a distance of R/D = 100 and angles ip. 



Figure 10. Comparison of predictions (Case SMC0008) with measurement of the SMCOOO convergent nozzle 
operating at Mj = 1.24 and TTR = 3.20 at a distance of R/D = 100 and angles ip. 
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Figure 11. Comparison of predictions (Case SMC00011) with measurement of the SMCOOO convergent nozzle 
operating at Mj = 0.90 and TTR = 3.20 at a distance of R/D = 100 and angles ip. 



Figure 12. Comparison of predictions (Case SMC00012) with measurement of the SMCOOO convergent nozzle 
operating at Mj = 0.70 and TTR = 1.00 at a distance of R/D = 100 and angles ip. 
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Figure 13. Comparison of predictions (Case SMC00014) with measurement of the SMCOOO convergent nozzle 
operating at Mj = 0.50 and TTR = 1.00 at a distance of R/D = 100 and angles i/j. 



Figure 14. Comparison of predictions (Case SMC0142) with measurement of the SMC014 convergent-divergent 
nozzle, Md = 1.185, operating at Mj = 1.185 and TTR = 1.80 at a distance of R/D = 100 and angles tjj. 
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Figure 15. Comparison of predictions (Case SMC0161) with measurement of the SMC016 convergent-divergent 
nozzle, Md = 1.50, operating at Mj = 0.784 and TTR = 1.00 at a distance of R/D = 100 and angles 'll). 



Figure 16. Comparison of predictions (Case SMC0163) with measurement of the SMC016 convergent-divergent 
nozzle, Md = 1.50, operating at Mj = 1.294 and TTR = 1.00 at a distance of R/D = 100 and angles 
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Figure 17. Comparison of predictions (Case SMC0164) with measurement of the SMC016 convergent-divergent 
nozzle, Md = 1.50, operating at Mj = 1.50 and TTR = 1.00 at a distance of R/D = 100 and angles ip. 



Figure 18. Comparison of predictions (Case SMC0165) with measurement of the SMC016 convergent-divergent 
nozzle, Md = 1.50, operating at Mj = 1.734 and TTR = 1.00 at a distance of R/D = 100 and angles ip. 
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Figure 19. Comparison of predictions (Case SMC0182) with measurement of the SMC018 convergent-divergent 
nozzle, Md = 1.80, operating at Mj = 1.80 and TTR = 1.648 at a distance of R/D = 100 and angles ■0. 



Figure 20. Comparison of predictions (Case Axi5) with measurement of the Boeing convergent dual-stream 
nozzle operating at Mj p = 1.281 &; Mj s = 0.848 and TTR p = 2.87 &; TTR S = 1.00 at a distance of R/D = 100 and 
angles i/j. 
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